Two-photon approximation in the theory of the electron recombination in hydrogen. 
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A rigorous QED theory of the multiphoton decay of excited states in hydrogen atom is presented. The "two- 
photon" approximation is formulated which is limited by the one-photon and two-photon transitions including 
cascades transitions with two-photon links. This may be helpful for the strict description of the recombination 
process in hydrogen atom and, in principle, for the history of the hydrogen recombination in the early Universe. 
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CN ! I. INTRODUCTION 

The recent accurate astrophysical observations and measurements of the cosmic microwave background (CMB) tempreture 
and polarization anisotropy Q|], $2^ triggered a new interest to the theory of the two-photon processes in hydrogen in view of 
. the important role of these processes in the cosmological hydrogen recombination. The history of the hydrogen recombination 
in the early Universe is described in many reviews, for example yD. The bound-bound one -photon transitions from the upper 
levels to the lower ones did not permit the atoms to reach their ground states: each photon released in such a transition in one 
atom was immidiately absorbed by another one. In particular, the Lyman-alpha radiation 2p-ls, being reabsorbed, reemitted and 
again reabsorbed, did not allow the radiation to escape the interaction with the matter. As it was first established in J4|], ||5|] the 
two-photon 2s- Is radiative transition presents one of the main channels for the radiation escape from the interaction with matter. 
Hence, the recent properties of the CMB are essentially defined by the two-photon decay processes during the cosmological 
recombination epoch. 

In la], tM it was argued that the ns — > Is (n > 2) and nd — > Is two-photon transitions can also give a sizeable contribution 
to the process of the radiation escape from the interaction with the matter. Recently this problem was investigated thoroughly 
in the theoretical astrophysical studies in 0. There is a crucial difference between the decay of the ns (n > 2) or nd 
levels and the 2s decay level. This difference is due to the presence of the cascade transitions as the dominant decay channels 
in case of ns (n > 2) and nd levels. For the 2s level the cascade transitions are absent. Since the cascade photons can be 
O I effectively reabsorbed, the problem of separation of the "pure" two-photon contribution from the cascade contribution arises. 
'— 1 . An interference between the two decay channels also should be taken into account. This problem appears to be not at all trivial 
and requires an application of the methods of the Quantum Electrodynamocs (QED) for the bound electrons. 

Quantum Mechanical theory for the two-photon transitions was first developed by Goppert-Mayer ill (III and the first evaluation 
of the two-photon 2s —> \s+2 r ){El) decay rate in hydrogen was performed by Breit and Teller [11]. The accurate nonrelativistic 
\Q \ calculation for this transition rate was given in 11211 : fully relativistic calculations, valid also for the H-like Highly Charged Ions 
(HCI) with arbitrary Z (nuclear charge) values were performed in II 1 311 — H 1 311 . The most accurate recent calculation for this 
transition rate with the QED radiative corrections taken into account belongs to Jentschura [16]. As well as for the neutral 
• hydrogen, the cascade problem does not arise for the transition 2s — > Is + 2j(El) in the HCI with arbitrary Z values. 

The two-photon transitions were investigated theoretically and experimentally also in the few-electron and many-electron 
atoms and ions. In particular, the two-photon transition ls2s 1 S'o — > (Is) 2 1 Sq + 2j(El) transition rate for the neutral He atom 
was first evaluated in 11711 . This decay channel also does not contain cascade contribution. 
J> , The cascade problem first did arise in connection with the decay of the metastable 2 3 Pq level in He-like Uranium: 2 3 Pq —> 
l 1 S'o+7(Sl)+7(Ml). In this case there are two possible cascade transitions: 2 3 P Q -> 2 3 S 1 +-f(El) l 1 S +-f(El)+'y(Ml) 
and 2 3 P -t 2 3 P X + 7 (Ml) -> l 1 ^ + j(M 1) + "/(El). The corresponding decay rate was first evaluated by Drake IU9ll . Later 
Savukov and Johnson 1 20] performed similar calculation for a variety of He-like ions (50 < Z < 92). In 01, JH the "pure" 
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two-photon contribution was obtained by subtraction of a Lorentzian fir for the cascade contribution from the total two-photon 
decay frequency distribution. A rigorous QED approach for the evaluation of the two-photon decay probability in presence of 
cascades was developed in ll2lll on the basis of the Line Profile Approach (LPA) in QED, i.e. the QED theory of the spectral 
line profile (see Ir22ll ). The LPA consists of a standard evaluation of the decay probability as a transition probability to the 
lower levels. In the presence of cascades the integral over emitted photon frequency distribution becomes divergent due to the 
singular terms, corresponding to the cascade resonances. To avoid such a singularity, the resummation of an infinite series of the 
electron self-energy insertions into the electron propagator was performed in [21]. This resummation converts into a geometric 
progression and in this way the electron self-energy matrix element (and the level width as its imaginary part) enters the energy 
denominator and shifts the pole from the real axis into the complex energy plane, thus making the transition probability integral 
finite. With this approach F. Low [ 23] first derived the Lorentz profile from QED. In lfl9tl . Il20h the level widths in the energy 
denominators were also introduced, though as the empirical parameters. In luDj the ambiguity of the separation of the "pure" 
two-photon decay and cascades was first revealed for HCI; it was shown also that the interference terms can essentially contribute 
to the total decay probability. 
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Nearly at the same time when the pap er lfl9ll did arrive, the cascade problem was discussed also for the ns (n > 2), nd 
transitions in the hydrogen atom ll24ll . ll25ll . In these works the "pure" two-photon contribution was obtained simply by omitting 
the resonant (singular) terms, responsible for the cascades. This approach was criticized later in |8]. Another ("alternative") 
method which formally allows for the separate determination of the "pure" two-photon contribution in case of the two-photon 
transitions with cascades was developed in a series of works by U. Jentschura 12611 - M29I1 . This app roach contradicts to the LPA 
results. The LPA was applied to the 3s — Is transition (including cascade) in hydrogen in [30], where the ambiguity of the 
separation of the "pure" two-photon and the cascade contributions was again demonstrated numerically, as in the case of the 
HCI j2lll . Very recently a paper lf3~lTl did arrive where an attempt was made to find a compromise between LPA and "alternative" 
approach. A reasonable agreement between the numerical results obtained by both methods was found. However, to our mind, 
the disagreement between the LPA and "alternative" approach is of conceptual character and cannot be eliminated. 

Thus from the QED point of view only the total two-photon frequency distribution has a direct physical sense in case of 
the two-photon decays with cascades. This quantity should be a basic tool for the description of the two-photon processes in 
astrophysics. The employment of the "1+1" approximation for the description of cascades should be avoided. Along this way the 
most recent astrophysical theories IH, are built. Still the "1+1" approximation is not fully excluded from the considerations 
in A, 0. 

In view of the recent very accurate (with relative accuracy ~ 1%) measurements of the properties of CMB |Q]], Jill and with 
expectation of the even more accurate (~ 0.1%) measurements in the near future, the theory of the cosmological recombination 
free of any uncertainties connected with the separation of the "pure" two-photon and cascade contributions should be formulated. 

In the present paper we will formulate such a theory for the two-photon and the multiphoton decays in hydrogen. In this theory 
only two types of the level decays should be present: the direct one-photon decays when they are allowed and the total two- 
photon decays without separation of the "pure" two-photon decays and cascades. The total solution of the problem formulated 
above consists of two steps. First, the pure QED problem of the description of the multiphoton transitions in hydrogen in the 
"two-photon" approximation should be resolved. That is, all the decays of the excited levels should be classified and described 
either as the direct one-photon transitions to the ground state or as the two-photon transitions with cascades. In the "two-photon" 
approximation transitions with more than two nonresonant photons should be neglected. The formulation of the "two-photon" 
approximation should finalize the first step of the studies. The present paper will concern only this first step. 

An important feature of the rigorous QED treatment of the process of recombination is that we have to trace the decay of 
every particular level up to the ground state. This is of course not the full picture of the recombination process. To be more 
close to the cosmological recombination one has to consider the transitions from the continuous states (plasma electrons) down 
to the ground state, taking into account the rescattering processes. This would correspond to the second step mentioned above. 
However, the accurate treatment of the recombination process from the particular excited level, as presented in this paper, also 
may be of interest. In particular, we demonstrate that the consequent QED treatment of the 3 p le vel decay should include the 
two-photon contribution comparable with the widely discussed two-photon decay of 3s level ll24ll - ll3lll . In this paper we limit 
ourselves only with electric dipole transitions (both in one- and two-photon decays) and ignore n'd — > ns transitions which also 
are of importance ([6]-[9[]). 

At the second step one should modify the basic astrophysical equations describing the level population in hydrogen in such a 
way that the imput data for them should be, apart from the direct one-photon transition probabilities, only the total two-photon 
decay rates, including cascades, without separating out the "pure" two-photon decay rates. The use of the "1+1" approximation 
should be fully avoided. This task is beyond the scope of our paper. 

Our paper is organized as follows. In Section II we formulate the basic concepts for the LPA-based theory for the two-photon 
decay with cascades. The two-photon approximation for the description of the multiphoton transitions in hydrogen is introduced. 
In Section III the standard derivation of the transition rate for the Lyman-a 2p — Is transition is presented and in Section IV 
the standard QED derivation of the Lorentz profile for this emission process is given. The same is done in Section V for the 
two-photon decay 2s — Is. The decay of the 3s level is considered in Section VI and the ambiguity of separation of the "pure" 
two-photon and cascade contributions is demonstrated. The decay of 3p level in the two-photon approximation is described 
in Section VII, where it is shown that this decay also contains the two-photon contribution comparable with the two-photon 
contribution to the decay of 3s-level. In Section VIII an investigation of the decay of 4s level in the "two-photon" approximation 
is performed which gives the clue to the general formulation of the two-photon approximation in the theory of the multiphoton 
transitions. Section IX contains discussion of the results and conclusions. 



II. TWO-PHOTON APPROXIMATION FOR THE MULTIPHOTON DECAYS WITH CASCADES. 

In this Section we will follow the derivation in [ 30] using this example for the formulation of the general principles of the "two- 
photon approximation" in the QED theory of the level decays. The grounds of this theory consist of few basic principles. First, 
all the decays should be traced up to the ground (stable) state. Within the "two-photon approximation" only such transitions can 
be defined unambigously. The two-photon approximation assumes that we take into account either direct (allowed) one-photon 
transitions from the excited level to the ground state, or the total two-photon transitions which end up also at the ground level. 
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The two-photon transition rates consist of the inseparably mixed "pure" two-photon transitions and cascades. The contribution 
of cascades is dominating and can be of the order of the direct one-photon transitions. The "pure" two-photon contributions and 
the corresponding interference terms between the "pure" two-photon contributions and cascades define the level of accuracy of 
our theory: two-photon approximation. Thus, we neglect the "pure" three-photon, "pure" four-photon etc contributions, taking 
into account, however, the cascade parts of the multiphoton decays. 

As the direct one-photon transitions as well as the cascade transitions we will consider only the allowed El transition dominant 
in the nonrelativistic theory. The order of magnitude of the corresponding transition rates parametrically equals to (in relativistic 
units) Wgi = Cg± ] ma{aZ) A , where m is the electron mass, a is the fine structure constant, Z is the charge of nucleus and 
C^j 7 ' is the numerical coefficient. In particular, for the Lyman-alpha 2p~ ls transition in H-like ions C^ ) (2p-ls) = 732.722. 

Having in mind the astrophysical applications of our "two-photon" approximation we have to classify different decay channels 
with respect to their contribution to the radiation escape from the interaction with the matter. One of the main channels for this 
escape was already mentioned above: this is the two-photon decay of the 2s state. The order of magnitude of the "pure" two- 
photon decay rate (in case of 2s-level, when the cascade contribution is absent, the "pure" two-photon decay rate coincides with 

the total one) is W^xlx = C E\E\ mc? ( aZ f & In case of 2s ->■ Is + 2 7 transition Cj^i O - ls ) = 24.7547. We should 
stress that in the present paper we do not investigate in detail the process of the occupation of the 2s state. The situations when 
metastable state appears as an intermediate state in the cascade processes of transitions from the upper levels, will be included 
in the treatment of the cascades. In particular, in our treatment in this paper the total decay rate of the 3p level, apart from the 
one-photon decay rate (3p — ls), incorporates also the two-photon decay rate H^ 7 jj (2s — ls). This happens due to the 

existence of the cascade transition 3p — > 2s + 7 — > ls + 37. Then r tot (3p) = + const • W^ 2 ^ (*) (see Section VI for 

details). In the standard treatment the transition rate of the upper link of this cascade (3p — 2s) contributes directly to the 

r*ot(3p): Ttot(3p) — Wjsi + Wgi\3p — 2s) (**). The difference between these two situations can be explained by the two 
different types of the experiment. The equality (**) describes the laboratory experiment, when the photon with the frequency 
uj(3p — 2s) is registered. Another situation occurs in the astrophysical context when it is important how fast an atom in an excited 
state will reach the ground level, i.e. the recombination will be accomplished. For a given cascade this depends on the slowest 
link of the cascade. In particular, for the 3p level the slowest link of the cascade 3p —> 2s + 17—5- ls + 37 is the lower link: 
2s — > ls + 27 and the total width T tot (3p) is defined by the equality (*). In the other words, one has to distinguish between the 
lifetime of a certain level A (equation (**)) and the lifetime of an excited state of an atom, provided that initially this atom was 
in the state A (equation (*)). 

We remind that we consider here the decay processes in one single atom and ignore the possibility of reabsortion of the photon 
with the frequency u) — E(3p) — E(2s) by another atom. 



HI. DECAY RATE FOR THE 2P LEVEL IN HYDROGEN. 



In Fig. 1 the decay scheme for 2p level (Lyman-alpha transition) is depicted. The emission process in frames of QED is 
described by the Feynman graph Fig. 2. The corresponding matrix element of the S'-matrix is given by (see, for example [32]) 



(A'\S W \A) 



d 4 x ip A > {x)~i^A* (x) ipA (x) 



(1) 



Here is the first-order 5-matrix, e is the electron charge, tp A (x) = VM(j)e %Ea *, tpA(r) is the solution of the Dirac 
equation for the atomic electron, Ea is the Dirac energy, 'ipA' = V'a'To is the Dirac conjugated wave function with %/A, being 
its Hermitian conjugate, 7^ = (70, 7) are the Dirac matrices and x = (f, t) is the coordinate 4-vector (r, t are the space- and 
time-coordinates). The photon field, or the photon wave function A fJi (x) looks like 



(2) 



where is the photon polarization 4-vector, k = (k, uj) is the photon momentum 4-vector (k is the wave vector, u) = \k\ is 
the photon frequency). 

The transition amplitude Ua'A is defined as 



(A'|5 (1) |A> = -2iri8(w - E A + E A ,)U$ A . 
Transition probability per time unit (transition rate) is defined via Ua'A as J3 



(3) 



W A , A = 2tt 



U 



(i) 



A' A 



5(lo-E a +E a >) • 



(4) 
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If the final state belongs to the continuous spectrum (as in our case due to the emitted photon) the differential transition 
probability should be introduced: 



dW A 'A(k,e) =2tt 



U (1) 

U A'A 



5(u-E A + E A , 



dk 
(2^)3 



(5) 



where dk = d 3 k = u! 2 dvdu>, dv is the element of the solid angle in the momentum space. Integration in Eq. (O over oj gives 
the probability of the photon emission with polarization e in the direction v = k/uj per time unit within solid angle dv: 



dW A > 



2n 



w A ' A 



-ikr 



A ' A 



dv . 



(6) 



where u) A i A — E A — E A > . The total transition probability follows from Eq. © after integration over angles and summation 
over the polarizations 



W A ' A = ^-u} A ' A ^ J dv [(e*a)e 



—ikr 



A ' A 



(7) 



For the atomic electron the characteristic scales for \r\ and |fc| = lj are: \r\ ~ 1/maZ, lu = E A > — E A ~ m(aZ) 2 . Then 
in the nonrelativistic case, in particular for the hydrogen atom (Z = 1), kr ~ a and the exponential function in the matrix 
element in Eq. can be replaced by 1. In the nonrelativistic limit the matrix element involving the Dirac matrices a (electron 
velocity operator in the relativistic theory) can be substituted by the matrix element of the operator p/m; where p is the electron 
momentum operator, with the Schrodinger wave functions. Then Eq. (0 takes the form 



w A , 



dv\(ep) 



A ' A \ 



(8) 



where the notation {...) A i A now implies evaluation of the matrix element with Schrodinger wave functions. Performing summa- 
tion over the polarization with the help of the standard formulas [32] and integrating over v yields 



4 e 2 

WjL'a = o— ua>a\(p)a>a\ ■ 
6 m 



(9) 



This is transition probability in the nonrelativistic limit in the "velocity" form. The "length" form W l A , A involving the electric 
dipole moment operator d = er of the electron can be obtained from Eq. (|9]l via the quantum mechanical relation 



u A ' A (r) A >A = —{p) A ' A 
m 



(10) 



Then 



W l A , 



A' A 



J A'A 



(d) 



A'A 



(ID 



Thus in the nonrelativistic limit only the electric dipole (El) photon emission is allowed. Using the atomic characteristic scales 
for \r\ and w, given above we easily obtain the order-of-magnitude estimates for the one-photon El transitions W 7 ^ 7 ' mentioned 
in the previous Section. In particular for 2p — Is transition (A = 2p, A' = Is) it follows from Eq. (fTTT i 



W^ ] (2p~ Is) = 732.722ma(aZ) 4 = 0.626 • 10 9 s" 1 



(12) 



IV. QED DERIVATION OF THE LORENTZ PROFILE. 



In this Section we give the QED description of the line profile. The basic ideas of the modern QED line profile theory were 
formulated by Low [23]. This method can be applied also for the description of the cascades and two-photon transitions. We 
employ the relativistic units (fi = c = 1) throughout this section. 
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A. Resonant scattering of the photon on the atomic electron. 

Consider first the elastic photon scattering process. The Fyenman graph corresponding to the process is presented in Fig. 3. 
According to the standard correspondence rules for the bound-electron QED (see, for example, HH) the S-matrix element, 
corresponding to the graph Fig. 3, is 

(A'\S^\A) = e 2 J d A Xl d A x 2 (^ A/ (x 1 h f , 1 A; i (x 1 )S(x 1 x 2 )^ 2 A; 2 (x 2 )^ A (x 2 )) , (13) 

where S(x\x 2 ) is the Feynman propagator for the atomic electron. In the Furry picture the eigenmode decomposition for this 
propagator reads (e.g. 13311 ) 

= ^J ^M) £ rf^tf > (14) 

2m J E n (l — iQ) + oo\ 

— oo 71 

where the summation in Eq. (fl4l i extends over the entire Dirac spectrum of electron states n in the field of the nucleus. 

Inserting the expression for the propagator in the S-matrix element, integrating over the time variables and frequency variable 
uj\ and using the connection between the 5-matrix and the amplitude U a>a Eq. ©, we obtain an expression for the scattering 
amplitude 

sc ^ E n -E A -u K " } 

with condition to = io' which implies the energy conservation. Here U u = ej^A^x) and lu denotes the frequency of the photon. 
In the resonance approximation the photon frequency u> is close to the difference of the two atomic levels: u> = u> res « Ea> —Ea- 
Then we can retain only one term in the sum over n in Eq. ( fT~5b : 

jj(2)res _ (Ulj)aA' A' A C lf c\ 

SC T71 771 

E A > - E A - L) 



B. Line profile for the emission process. 

It follows from Eq. ( fToT l that in the resonance approximation the emission amplitude can be expressed like 

U - (^)aa' n7 -, 

E A > - E A - to 

The absorbtion amplitude can be presented similarly. Expression ( TTol ) for the scattering amplitude is singular at the resonant 
frequency. To remove this singularity one has to consider the electron self-energy insertion in the electron propagator in Fig. 
3. According to ll23ll this implies the arrival of the Lorentz line profile for the scattering process. The lowest-order electron 
self-energy insertion is shown in Fig. 4. Using the correspondence rules, we obtain an expression for the correction to the 
scattering amplitude 

^{K) Ani \t{E A + u)] (U u ,) nsA 

jjii) = _ iliilis (18) 

{E ni -E A -u){E n3 -E A -u) ' 1 ; 



where is the electron self-energy operator for the bound electron [32 

In the resonance approximation n\ = n 2 = A' and the correction to the scattering amplitude is 

7-/(4) _ _TT{2)res\ J A' A' /iqn 

Ea'-Ea-u ' 

Repeating the insertions in the resonance approximation we obtain the geometric progression. The summation of this progression 
yields 

jjres _ (Uu)aA' (Uu')a'A 

E A > - E a -uj 
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where E A > = E A + (t{E_ 
expression 



a + ) . The emission amplitude in the resonance approximation is then presented by an 

A' A' 



uiJAA' 



E A > - E A - lo 



(21) 



Apart from the electron self-energy (SE) also the vacuum polarization (VP) insertion in the electron propagator in Fig. 3 should 
be considered to all orders in the resonance approximation. The lowest-order VP insertion is described in Fig. 5. The VP 
insertions lead to the following change of the nenergy denominator: 



E A > =E a +(£(E a +lu 



A ' A ' 



n 



A ' A' 



(22) 



where (II) A i A i is the vacumm polarization operator for the bound electron 13211 . 



The real part of the matrix element (y>(E/ 



i! A + w) ) presents the lowest order contribution of the electron self-energy to 

/ A' A' 

the Lanb shift, the imaginary part of this matrix element defines the radiative width T A i of the level A': 



Z(E a +lu) 



— Li ai 



SE 



'A.' A' " 2 

The other lowest-order part of the Lamb shift is the vacuum polarization part 

(TL)a'A> = 

The vacuum polarization does not contribute to the width. Then 



r A > 



tVP 
L A , . 



(23) 



(24) 



E A . = E A . + L A , - -Y A , , 



where L A > — L°Jr' + L V A T and the emission amplitude looks like 



E A , +L A -E A -u- %r A , 



The total transition probability A' — > A is 



dW AA > (u) = — ^ J dv\U ern \ 2 uJ 2 doj . 



Insertion of the expression d26b in Eq. ( |27] | in the resonance approximation yields 

i r AA >dui 



dW AA , (u) 



2tt (E a , + L A , -E A - u>y + \T\, 



(25) 



(26) 



(27) 



(28) 



Here T AA i is the partial width of the level A' connected with the transition A' — > A. Equation d28l) defines the usual Lorentz 
profile for the emission spectral line. The resonance frequency in zero-order approximation is defined by lo — uj res w E A i — E A . 



V. TWO-PHOTON DECAY RATE FOR THE 2s Is TRANSITION. 



In this section we discuss the two-photon processes 2s — > Is + 2j(El) which corresponds to the decay scheme Fig. 6. The 
two-photon transition probability A — > A' + 2j corresponds to the second-order 5-matrix element Eq. ( fT3l > (see Fig. 7) with 
the replacement of the absorbed photon by another emitted one. Using again Eqs. ([3]) and for the two-photon transition and 
integrating over time and frequency variables in Eq. ( fT3b . we find for the sum of the contributions of the both Feynamn graphs 
(see Figs. 7) 



dW A 'A = 2tt(5 (E a -Ea'-w- oj') 



U A'A 



dk dk 1 
(2tt) 3 (2tt) 3 ' 



(29) 
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where 



(2) _ 2ire A 

U A'A — I 



(cLAlS) (al% p ) ^Uaz,-.,) (all-.) 

\ e >«/ A'n V e > k / nA , V e > k / A'n V e > k J nJ 



^ E„ — Ea + U>' ^ E n — E t 



(30) 



A g £ = ee lkr and a are the Dirac matrices.. 

In what follows, we will be interested in the decay rate of the ns levels (A = ns, A' = Is) in hydrogen. In this section we 
focus on the case n — 2, when the cascades are absent. The schematic picture of the decay 2s — s- Is + 2j(El) is given in Fig. 
6. In the nonrelativistic limit, after the integration over frequency uj', over photon directions dv, dv 1 and summation over all 
polarizations e, e ', we obtain for the photon frequency distribution: 

g^j3/^^ Uj)^ 2 

dW2s,x s (uj) = — e 4 |Si Sl 2 S (w) + S , i Si2s (^o - w)| 2 dui, (31) 

Z ITT 

g (u) = y (Ris\r\Rn' P )(Rn> P \r\R2s) (32) 
, E n i p — E ns + uj 



(Rn>i>\r\Rni) = J r 3 R n , v {r)R nl {r)dr , (33) 



where wo = E2 S — Ei s , R n i(r) are the radial part of the nonrelativistic hydrogen wave functions, and E n i are the hydrogen 
electron energies. Here we have used again the quantum-mechanical relation Eq. ( fTOb ; Eq. OTb is written in the "length" form. 
The decay rate for the two-photon transition can be obtained by integration of Eq. ( f3TT > over the entire frequency interval 

W 2s ,i s = ~ J dW 2s ,i s (u)). (34) 
o 

In case of 2s state the cascade transitions are absent, the frequency distribution Eq. (l3TT l is not singular and the integral Eq. (l34T > 
is convergent. The result of the integration over frequency ui is well known and equal to W^s.is = 24.788ma 2 (aZ) 6 r.u. = 
8.229 s" 1 lfl2t 

VI. TWO-PHOTON DECAY WITH CASCADES FOR THE 3s - Is TRANSITION. 

In case of the cascade transitions (n > 2), some terms in Eq. ( l32t become singular and the integral Eq. d34l l diverges. This 
divergency has a physical origin: an emitted photon meets the resonance. The corresponding scheme of the decay for n = 3 is 
given in Fig. 8. So the divergency can be avoided only by introducing the width of this resonance. This situation was studied 
in ll2lll for the HCI. The same recipe can also be used in case of the hydrogen atom. Following the prescriptions given in 1I2II1 
we separate out the resonant terms (corresponding to cascades) in the sum over the intermediate states Eq. (f32l > and apply 
Low's procedure i23ll for the regularization of the corresponding expressions in the vicinity of the resonance frequency values. 
Practically this leads to the apperance of the energy level widths in the energy denominators. Then the Lorentz profiles arise for 
the resonant terms in the expression for the probability. However, the Lorentz profile is valid only in the vicinity of the resonance 
and cannot be extended too far off from the resonance frequency value. As for any multichannel processes such a separation is 
an approximate procedure due to existence of the interference terms. 

The integration over the entire frequency interval [0, ujq] in Eq. (|34l should be split into several subintervals, e.g. 5 in case of 
the two-photon emission profile for the 3s-level decay, see Fig. 9 yOfl . The first interval (I) extends from uj = up to the lower 
boundary of the second interval (II). The latter one encloses the resonance frequency value lu\ = E^s — Within the interval 
(II) the resonant term n = 2 in Eq. d32l should be subtracted from the sum over intermediate states and replaced by the term 
with modified energy denominator (see Section IV). This modified denominator is E^ v — Es s +lj + where T — T 2 p + ^3 S - 
The third interval (III) extends from the upper boundary of interval II up to the lower boundary of the interval (IV), the latter one 
enclosing another resonance frequency value CJ2 = E 2p — E\ s . Within the interval (IV) again the resonant term n = 2 in Eq. 
J32l should be replaced by the term with modified denominator E 2p — E\ s — uj — ^T2 P - Finally, a fifth interval (V) ranges from 
the upper boundary of the interval (IV) up to the maximum frequency value luq. Note, that the frequency distribution dWj, s _\ s (uj) 
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is symmetric with respect to w = ojo/2 with a 1% accuracy (the asymmetry is due to the difference between V — I^p + T 2s and 
r 2p , respectively). 

Inserting Eq. ( f3TT > into Eq. d34l i and retaining only the resonant term within the second and fourth frequency intervals, will 
yield the cascade contribution to the total two-photon decay rate of the 3s-level. Taking the ratio to the total width of the 3s-level 
T 3s we will obtain the absolute probability or branching ratio W^ a ^ a /y 3s = fr^^p-is for the cascade transition. The 

contributions to o^^p-ls fr° m me intervals (I), (III), (V) are assumed to be zero. The cascade contribution of the 3s-level 
results (in the "length" form) 



jrr(cascade I7) 
" 3s;ls 



(ii) 



{R 3s (r)\r\R 2p (r))(R 2p (r')\r'\R ls (r')) 



— e 3s + u — §r 



doj + 



(35) 



-— J U 3 (U Q -L0)3 
(IV) 



[R 3s (r)\r\R 2p (r))(R 2p (r')\r'\R ls (r')) 



E2p — Ei s — <*> — 2~^2p 



du. 



According to the discussion in Section V the "pure" two-photon decay probabilities within each interval, defined in Section 
V, look like 



dw^; 2j) - ^u\u Q -uf - 



du, w € II 



(36) 



dW^T 7) = =-u\u -uf Si.i3.M + 5i^,(a;o-a;) 



27tt 



j(2p) 



du, u € IV 



(37) 



27-k 



u 3 (u - u) 3 \Si s . 3s {uj) + Si s - 3s {u - u)\ du, weI,III,V. 



(38) 



Here S[ 2 s p 3s (u) is the expression (l32b with the n — 2 term being excluded. 

Unlike cascade, all the intervals contribute to the "pure" two-photon transition. The branching ratio for this transition 3s 
27 + Is appears to be 



, (puro27) 
J 3s-ls 



1 1 
2~T~ 



3s 



dW. 



(pure2"f) 



3s:ls 



( W ). 



(39) 



It remains to introduce the interference contribution. This contribution comes only from the intervals II and IV. The correspond- 
ing frequency distribution functions are given by 



(inter) _ 4a; 3 (w - u)' c 



27n 



Re 



(R 3s (r)\r\R2 P (2r))(R2 P (r')\r'\R ls (r')) 

E 2 p — E 3s + UJ — ip 2p 



du , u e II (40) 



dW. 



(inter) _ 4o; 3 (a;o - uf 



3s;ls 



27tt 



Re 



(R 3s (r)\r\R2p(2r))(R2p(r')\r'\R ls (r')) 



2l> 



2 1 2p 



S ls]3s (u) + s[ 2 s p 3 ) s (u - u) 



du, ue IV (41) 



and branching ratio results as 



, (inter) _ J- 

3s;ls 2p 



3.s 



dW- 



(inter) 1 



3s;ls 



2T 



•3 s 



dW. 



(intor)2 



3s;ls 



(42) 



(ii) 



(IV) 



The results of our calculations are presented in Table 1. It is convenient to define the size Aw of the second interval as 
multiples I of the widths T, i.e. Au = 21T and for the fourth interval as Aoj = 21T 2p , respectively. In Table 1 numbers are given 
for different values of I ranging from I ~ 10 5 up to I ~ 10 7 . The upper boundary of interval II equals U\ + IT = + IT (in 
a.u.), while the lower boundary of interval IV equals 0J2 — lT 2p = § — T 2p . The different lines of the Table 1 present branching 
ratios and transition rates of the "pure" two-photon and "interference" channels, respectively. For the more detailed analysis the 
contributions of the "pure" two-photon transition rate for the each frequency interval are also compiled. The branching ratio 
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and the transition rate for the cascade contribution can be obtained from the relation 63^20— i« + ^2s-is 7 + ^Is-is = 
relation is sutisfied with high accuracy since the only decay channel neglected is the very weak direct 1 -photon Ml transition 
3s — > Is + 7. From the Table 1 we can draw the same conlusions: as in the case of the HCI 1I2TII : the "pure" two-photon and 
cascade contributions to the total decay rate appear to be inseparable. Changing the interval size Aw, we obtain quite different 

values for d^f 7) ranging from 202.16 s" 1 (for I = 10 4 ) up to 7.9385 s" 1 (for I = 1.00256 • 10 7 ). 

Moreover, in our calculations - depending on the size of the interval - the interference contribution also can become quite large, 
comparable in magnitude with the "pure" two-photon contribution. Thus, we demonstrated that even the order of magnitude of 
the "pure" two-photon decay rate for the 3s-state in hydrogen can not be predicted reliably. 

Earlier the result 8.2196 s _1 for the "pure" two-photon decay of the 3s-level was reported in l24ll and confirmed in l25ll . 
However, as it was pointed out in H] in both papers f24ll . l25ll the summation over the intermediate states was not performed 
properly. The "nonresonant" contribution 10.556 s" 1 deduced in @], which plays the role of the "pure" two-photon decay rate 
is well within the range of our values given Table 1. However, the result 2.08 s _1 obtained for the "pure" two-photon decay rate 
in lElll is in strong contradiction with the present analysis (see the discussion in ll30lo . 

Very recently, a paper lf3~lll did arrive where both the standard QED approach, based on the line profile theory ( lfl9ll - ll2~lll ) and 
the "alternative" approach based on the two-loop Lamb shift theory ([26]-[29]) were applied to the calculation of the two-photon 
transition in hydrogen. A reasonable agreement between the two methods was found. However, from the derivations in 113011 it 
follows that the employment of the Lamb shift imaginary part gives exactly the same results as the LPA QED approach. 



TABLE I: Branching ratios and transition rates (in s x ) for the different decay channels for the decay probability of the 3s level with different 
frequency interval size (/)• 



I 


10 4 


10 5 


2.5 ■ 10 5 


5 ■ 10 5 


10 6 


1.5 ■ 10 6 


4.53 ■ 10 6 


1.00256 ■ 10 7 


^(pure27) 


3.2003 • 10" 5 


3.5091 ■ 10~ 6 


1.6270 • 10 -6 


1.0239 ■ 10~ 6 


7.6765 • 10" 7 


7.2201 ■ 10~ 7 


9.1487 • 10 -6 


1.2567- 10~ 6 


py(purc27) 


53.054 


7.0547 


3.5743 


2.1898 


1.27737 


0.85130 


2.4979 ■ 10~ 6 





^(pure27) 


0.006247 


0.06247 


0.15614 


0.31201 


0.62183 


0.92718 


2.4666 


3.9810 


TJ/(pur c2 7) 
VV III 


95.536 


7.8778 


2.7928 


1.4517 


1.0457 


1.0031 


0.86005 





TJ/(P urc2 7) 
VV IV 


0.006185 


0.061847 


0.15458 


0.30890 


0.61569 


0.91813 


2.4523 


3.9575 


jy(purc27) 


53.561 


7.1101 


3.5999 


2.2056 


1.2886 


0.861254 


3.1665 ■ 10 -4 





T^/-(purc27) 


202.16 


22.167 


10.278 


6.4680 


4.8492 


4.5609 


5.7792 


7.9385 


jointer) 


-1.4342 ■ 10~ 9 


-1.4343 ■ 10" 8 


-3.5852 ■ 10" 8 


-7.1665 ■ 10~ 8 


-1.4302 ■ 10" 7 


-2.1376 ■ 10~ 7 


-6.0829 ■ 10~ 7 


-1.0459 ■ 10~ 6 


^(intcr) 


-0.0090599 


-0.090602 


-0.22647 


-0.45270 


-0.90346 


-1.3503 


-3.8426 


-6.6067 



In the end of this Section we will explain why the contributions of the "pure" two-photon transition rates in Table 1 are of the 
same order as the interference terms. As it was mentioned earlier the contribution of the cascade is of the order ma(aZ) 4 in r.u. 
However, this is the result of the integration over the frequency interval of the order I^ 17 ) ~ ma(aZ) 4 . Then the order of the 
magnitude of the "amplitude" in the integrand is 1 (and dimensionless). The order of the magnitude of the "pure" two-photon 
contribution is ma 2 (aZ) 6 r.u. This result also incorporates the integral over frequency interval of the order u>o = m(aZ) 2 in r.u. 
(seeEq. (34)). Then the order of the magnitude of the corresponding "amplitude" in the integrand is [ma 2 (aZ) 6 /m(aZ) 2 ] 1 ^ 2 = 
a(aZ) 2 . The latter value is again dimensionless. Multiplying the cascade and "pure" two-photon "amplitudes" in the integrand 
we will have the dimensionless integrand of the order 1 • a(aZ) 2 . Finally, integrating the product over the frequency interval 
r ~ ma(aZ) i , where the interference terms are nonzero, we obtain the contribution of the order ma 2 (aZ) i r.u. This is 
parametrically the same as the order of magnitude for the "pure" two-photon contribution. 



VII. "TWO-PHOTON APPROXIMATION" FOR THE TREE-PHOTON 3P-1S TRANSITION IN HYDROGEN. 

The 3p — Is decay can occur either as one-photon or as three-photon process. These channels do not interfere due to the 
different number of photons in the final state. The one-photon decay 2>p —> Is + j(El) corresponds to the decay scheme Fig. 

1, where the initial state 2p should be replaced by 3p. The value (3p — Is) can be obtained from the formula (11) when 

inserting there A = 3p, A' = Is. The result is 

W^\3p- Is) = 195.613ma(aZ) 4 = 1.67342 • 10 8 s" 1 . (43) 

The schematic picture for the process A — » A' + 37(^1) is given in Fig. 10 and the corresponding Feynman graph is depicted 
in Fig. 11. The three-photon emission probability was evaluated in our work lf3~4h for the 2p — Is three-photon transition 
in hydrogen. The parametric estimate can be easily found and is equal to ma 3 (aZ) 8 in relativistic units. The value of the 
probability for this transition is 0.4946ma 3 (a!Z) 8 r.u. 
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The 2p —> Is + 3^f(El) transition is "pure" 3-photon transition (no cascade are possible in this case)and, according to our 
general scheme (see Section II) should be neglected. For 'dp — > Is + 37(£'l) transition the two cascades should be taken into 
account: 3p — s- 2s + j(El) Is + 3j(El) and 3p — >• 2p + 27 (.El) Is + 37(E1). The contribution of these two cascades 
will be studied in this Section. It will be shown that the decay rates of these two cascade channels will be comparable with the 
"pure" two-photon contribution of the 3s — > Is + 2j(El) channel. Unlike the 3s — > Is + 27(E1) case, where the contribution 
of the "pure" two-photon decay is nonseparable from the cascade contribution, for 3p — > Is + 37(E1) decay only the cascade 
contribution should be taken into account. The reason is that the "pure" 3-photon contribution is beyond the accuracy of the 
"two-photon" approximation, adopted in this paper. 

The S'-matrix element for the 3-photon decay process A —> A' + 3j is: 

S*A = H) 3 / d i x 1 d i x 2 d 4 x^ A ,{x 1 ) (7^?" (zi)) S{x lX2 ) (7^^(3:2)) S(x 2 x 3 ) (^ 3 A™(x 3 )) ^ A (x 3 ) , (44) 

where the photon field A*" (x) is described by Eq. d2j and ui, uj' uj" denote the frequencies of the photons. The electron 
propagator S{x\x 2 ) and the electron wave functions ip A > (%x), iPa(x3) are defined like in Section EH: 

$a>(xi) = i>A>{n)e iEA ' tl ; 4j a (x 3 ) = 4>A(?3)e- lEAt \ 

Using Eqs (|2), ( fl4l > and (|45T > we can perform the time integration over time variables in Eq. d45l > 

J dt 1 e i( - EA ' +u " +U1 ^ 1 = 2tt8 {E a > + u" + wi) , 

dt 2 e 1 ^ '-"i+^)ta = 2n5 (uj' -u 1 + oj 2 ) , 

dt3e i( U - U2 -E A )t 3 = ^ _ ^ _ 



(45) 



(46) 



Then the frequency variables in the two energy denominators are ui\ = —E A > — uj", ui 2 = u>x — uj' = —Ea — uj" — uj'. From 
these two equations follows lu + uj' + uj" = Ea — Ea', what is the energy conservation law for this process. 
Then after the integration over u>i and iu 2 the S'-matrix element can be written in the form: 



(A'\S^\A) = (-l) 3 e 3 



d 3 r 1 d 3 r 2 d 3 r 3 ijA'(r 1 )(e"a 1 )\j —e 

to 



-i(k"r 1 ) 



E 



ipm (n)ip ni (r 2 ) 



E ni (l-iO)-E A ,-u/< 



(47) 



(eWV^E 

V LJ ^— ' 



^2(^2)^2 (r 3 ) 



E n2 {l-iO)-E A 



7 (SSs)J— e~*^MWE A , - Ea+oj + uj' +uj"). 



Here Si (i = 1, 2, 3) is the Dirac matrix, k is the wave vector of the corresponded photon (\k\ = uj) and e is the polarization 
vector of the emitted photon. 

The transition probability per time unit is defined via Eqs (01, © by 



Wp2 = 27T 



u 



(3) 
A' A 



5(E A > - E a +uj + uj' + uj"). 



The differential transition probability is introduced like 



dW A'A (*> e; k', e'; k" ', e 7 ')2n 



rr(3) 
U A'A 



8{Ea> — Ea + uj + uj' + uj") 



d 3 k d 3 k' d 3 k" 
(2tt) 3 (2tt) 3 (2tt) 3 



where d 3 k = uj 2 duidv and dv is the element of the corresponding solid angle in the momentum space. 

The total transition probability follows from Eq. d49b after integration over angles and summation over polarizatoins 



W 



(37) 



A' A 



3!(2tt) 5 J^ 7i 



dv / dv' / dv" / uidui / u'duj' / u>"du>" x 



E 



(A'|(e 7 'a)e- l ( fc " ?1 Vi>(^il(^a) e "' (P?2 V2)(^2|(eS)e- l ^ f ' 3 )|A) 
(E ni - E A , - uj")(E n2 - E A , -uj'- uj") 



(48) 



(49) 



(50) 



where k, e; k', e'; k", e" are the wave vectors and the polarization vectors for the three emitted photons. In Eq. (|50l the 
permutation symmetry of the emitted photons is taken into account. 
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Consider now the 3p — > 37(£1) + Is transition. In this case the initial state A = 3p and final state is A' = Is. Due to the 
conservation law we can write uj' + uj" = E 3p — E ls — uj and, therefore, the transition probability (in dipole approximation) is 

W[l% = ^-y; j dP j ^ J dl? ' J ^ J " W " W") X 

(ls| (e"p) |m) (n x | (e'p)|n 2 ) (n 2 1 (ep) |3p) 



E 



ni ,ri2 



- - uj")(E n2 - E A , + u) 



(51) 



where ujq = E^ P — E\ s and p is the electron momentum operator. 

This three-photon emission probability contains resonant transition, when E n2 = E 2s . The value E„ 2 — E 2s corresponds to 
the pole in the integral over frequencies in Eq. (|5H uj = E-^ p — E2 S . Then wo — uj — uj" = E^ p — E\ s — E^ p + E 2a — uj" = 
E 2s — Eis — uj" at the pole value. The integral over uj' becomes divergent. It can be regularized in the resonance region according 
to Section IV. After the regularization procedure we can integrate Eq. dBTT l over frequency uj: 

v u mmmf . eS^. (52) 

The T^ p 2 2s i s tne partial width corresponded to the 3p — > I7 + 2s transition and T^p is the total width of the 3p level which is 
the sum of the all partial widths for the transitions from 3p state to the ground state. 

The remaining factor in Eq. ( f5Tb represents the two-photon 2s — Is transition probability. Then the expression for the resonant 
contribution after collecting all the terms in the right-hand side of Eq. dBTT l reads 

n r (l7) 
w (37)resl _ 3 M/ (2 7 ) 1 3 P -2s 

"3p-ls — ~7 yv 2s-ls~^ ' \ D3 ) 

4 1 3p-l s 

In the sum over all intermediate states (ni) in Eq. ( BTT l exists also another resonant term: E ni = E 2p , when ui" = E 2p — Ei s . 
In this case ujq — ui — uj" = E^ p — E\ s — E 2p + Ei s — uj = E^ p — E 2p — uj. Regularization in the region uj" = Ez p — E 2p 
leads to the final expression 

wr (3 7 )re«a 3 (2 7 ) 

Thus in the "two-photon" approximation the total probability of the 3p level decay can be expressed as 

W^-L = W$? u + \wg% + (55) 

The first term in Eq. d55l ) is the ordinary one-photon width. The two other terms represnt "cascade two-photon" contributions 
which are of the same order of magnitude as the "pure" two-photon contribution to the 3s level width T^ s . However, unlike the 
case of 3s level, where the "pure" two-photon contribution cannot be distinctly separated out ll30ll . the "cascade two-photon" 
contribution to T^p is given explicitly by the two last terms in the right-hand side of Eq. d55l ). The corresponding numerical 
values for hydrogen atom are 

o 

1 W al-2 P = 0.034005675 s" 1 

q W( 1 7) 

±^Ezl± W gZ\ s = 0.730334 s' 1 (56) 

4 r 3p 

which is not negligible in comparison with the order of magnitude of the two-photon contribution ~ 10 s _1 to the decay rate 
of 3s level. Evaluation of W^} 2 p m ^1- *E3 i s performed in the nonrelativistic limit, when the sublevels 2pi/ 2 and 2p 3 / 2 are 
degenerate. The calculation of the transition rates between separate fine-structure components of 3p, 2p levels is given in the 
Appendix A. The interference term between the second and the third terms in the right-hand side of Eq. d55l ) are absent as well 
as the interference with the one-photon transition (the first term in the right-hand side of Eq. d55l). The reason is that the second 
and third terms are nonzero close to the two different frequency values uj in the integral in Eq. (51). Therefore the product of the 
two corresponding amplitudes in the integrand in Eq. (51) is always small. 
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VIII. "TWO-PHOTON APPROXIMATION" FOR THE FOUR-PHOTON 4S-1S TRANSITION IN HYDROGEN. 



Consider now the transition 4s — > Is. Apart from the ordinary two-photon processes 4s — > 3p + 7 — > Is + 27 and 
4s — > 2p + 7 —> Is + 27 there are more complicated 4-photon processes with two-photon links. The schematic picture for the 
two-photon decays of 4s state is given in Fig. 12 and the 4-photon decay picture is presented in Fig. 13. In general, the S-matrix 
element for the 4-photon decay can be described by 

Sit]is = H) 4 / d i x 1 d 4 x 2 d i x 3 d' i x^ A ,(x 1 ) (tmi^T (a*)) S(x lX2 ) ( 7/12 ^f (xa)) * (57) 

S(x 2 x 3 ) (7^3 A %' (^3)) S^a^) (Tau^W) , 

where all the notations are the same as in Eq. (44). The electron wave functions ipA'(xi), ipA(xi) are the same as in Eq. (45), 
with A' = 4s, A = Is. Perfoming time integration and integration over frequencies w\, lj 2 , W3, we recieve 



(A'\S^\A) = (-l) 4 e 4 / d 3 r 1 d 3 r 2 d 3 r 3 d 3 r 4 i'A'(ri)(e T "d 1 ] 



-i(k"'r\ ) 



E 



ip ni {r 1 )ip ni {r 2 ) 
E ni (l-iQ)-E A r-u'- 



(^ 2 )W±L e -^) E 



V w' ^— ' 



£ n2 (l-iO)-^ 



^3 (fs^nsfo) 



£ n3 (l-iO)-^ 



7 x (58) 

1* 



-i(fcr*4) 



^A{r^)5{EA' — Ea + w + w' + w" + w'")< 



Here fc, e; fc', e'; fc", e"; fc'", e'" are the wave vectors and polarization vectors for the four emitted photons. Using Eq. (0) we 
can write the transition probability like 



wi 4 ]) = — 

AA' 4 , (2 ^ 



dv I dv' I dv 



dv" 



Ludu lo du uj dm I lj duj 



U A'A 



(59) 



Now we will demonstrate how the two-photon emission is included in this four-photon process. First, we consider the 
4s — >• 3p + 7 — > 2s + 27 —> Is + 47 transition. In this case we fix E n3 = E 3p and E ri2 — E 2s , then the resonant frequencies 
are ui = E± s — E$ p , 10' = E% p — E 2s . Therefore, in the dipole approximation, 



u ls4s 



E 



(ls\e"'p\n 1 )(n 1 \e 7 'p\2s)(2s\e'p\3p)(3p\ep\4s) 



(E ni - E 2s + u"){E 2a - E 3p + uj'){E, 



3p 



E is + u) 



other 23 terms 



(60) 



where 'other 23 terms' differ from the first one by permutation of matrix elements. We again use the expressions for the one- 
photon probabilities: 



(17) 



W, 



4s -3p 



(17) 
3p-2s 



^E/ dmp\ms) 

e 

^E/ dv'\{2s\e>p\3p) 



(61) 
(62) 



and the equality 



/ 



did 



(UJ + + ±1* 



2 arctan ^ 
f 



7T 

r 



(63) 



which is valid for the small (compared to ujq) values of T. 

Then, we can get for the four-photon transition probability an expression: 



W, 



(47) 



4s — Is 



4!(2tt) 



E 



e"e"' 
dv" I dv 1 ' 



(E 2s -E 3p +ojr + k T l 



w£l\ p <^ 



{E 3p -E is +uy + \Tl p 



E (ls|e>> 1 )(n 1 |^l2s) thcr23terms 
^ E ni - E 2s + ui" 



(64) 
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or, in equivalent form, 



W 



4.-1- 4!(2^)5 ^ y (^.-£3p + w') a + Jri y (E 3p -E 4s +oj) 2 + iT 2 3p 

e"e"' 

ls\e r "p\n 1 )(n 1 \e 7 'p\2s) (ls|e 7 'p|ni) (m |e"'p|2s) ' 2 



12 2 / <fz/' / di/ 



u> au) I us 



E 



- ^2 S + w" ^ E ni - E 2s + uj* 

Til tlx 



(65) 



and 



Finally, for the 4s — > 3p + 7 — > 2s + 27 — > Is + 47 emission process we get 

T „(4 7 )rcsl _ J ^4£-3p vv 3p-2s „ A (2 7 ) 

KK 4s-ls — o F ^s-ls ■ l. ot V 

1 4s 1 3p 

Also the processes 4s — > 3s + 27 — > 2p + 3j — > Is + 47 and 4s — > 3p + 7 — >■ 2p + 37 — )■ Is + 47 should be considered. 
Performing similar calculations we get 

w (4 7 )rcs2 _ " W 3s-2p W 2p-ls w (2 7 ) _ d ^3s-2p w (2 7 ) 

vv is-ls — o w 4s-3s = o ^ rK 4s-3s \° 1 ) 

£ 1 3s 1 2p ^ 1 3s 

q M/( l7 ) W (l7) "? T/l/ (l7) 

wi^r - 1 r 3 p ; p - u w&\ = \^^wg% p . (68) 

* 1 4s 1 2p Z 1 4 S 

Finally, for the four-photon decay of the 4s hydrogenic state within the two-photon approximation we recieve 

Wtotal _Tl/(27) 3 ^ 3 s-2p Ty (2 7 ) 3 ^4s-3p „,(2 7 ) 3 ^ 4 s-3p ^3p-2s „,(2 7 ) 

^45-15 — ^4s-ls + ^4s-3s + W 3p-2p + 9^ p ^s-ls • ( W > 

* 1 3s Z 1 4 S Z 1 4 S 1 3 p 

where wjfZis is the "pure" two-photon contribution to the 4s — ► Is decay (see Fig. 12). We remind that, in principle, this 
contribution is inseparable from the cascade contribution in Fig. 12. The numerical values for the cascade contributions with the 
two-photon links (three last terms in the right-hand side of Eq. (69)) are: 



q W^f) 

^^s-2p (2 7) ^ 0.00438791s" 1 



.} p ■ ■ 4s— 3s 
Z 1 3s 

(17) 



3 W^-3p w gy) = 0.0284699 s" 1 ( 70 ) 

2 T4s 

| ^^ ^1= 0.611441s- 1 

* 1 4s 1 3p 

These values are again not neglidgible compared to the contribution W^Z\ S . The latter can be estimated as ~ 12 s -1 (H]. 

For the evaluation of the 4s — > 3s+27 transition probability it is enough to consider the standard expression for the distribution 
function dW(u>) (Eqs d3T1l-(l34l) with the replacement of the wave functions 2s — > 4s, Is — > 3s. Note, that there are no 
intermediate p-states between 4s and 3s levels. The existence of the degenerate 4s, 4p and 3s, 3p levels does not change the 
result. The result of the calculations in the nonrelativistic limit is 

w£l\ s = 0.00879957a 2 (aZ) 6 r.u. = 0.00292527 s- 1 . (71) 

The relativistic value is equal 0.002924794 s _1 . The relative difference of relativistic and nonrelativistic values is about 0.016%. 
Similarly the nonrelativistic result for the 3p — 2p two-photon transition rate can be obtained: 

W 3p-2 P = 0.13618a 2 (aZ) 6 a.u. = 0.045271s" 1 . (72) 

For the summation over the nonrelativistic hydrogen spectrum in the expressions for the two-photon transition rates we used 
the Green function method (see, for example, JUt], l36ll and also |H3]). Here we give also the numerical result W^gl\ a = 
0.194113q; 2 (q;^) 6 r.u. = 0.0645296 s _1 with relative difference with the relativistic one 0.0082%, though this value does not 
enter directly to the 3p, 4s decays. 
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IX. CONCLUSIONS 

In our paper we have considered the processes of multiphoton transitions for hydrogenic atom. Recent astrophysical inves- 
tigations necessitate the detailed analysis of the multiphoton emission processes and, namely, of the separation of the "pure" 
multiphoton radiation and cascade emission. The "pure" two-photon emission leads to the photon escape from the matter and, 
thus, presents the formation mechanism for the background radiation. 

We began with the standard QED derivation of the Lorentz profile for emission process. After this we have investigated 
the decay of the 3s level and showed the ambiguity of separation of the "pure" two-photon and cascade contributions. We 
demonstrated that the strict separation of the "pure" two-photon and cascade contributions for 3s-level decay in hydrogen is 
impossible. Moreover, we show that even the approximate separation of these two decay channels cannot be achieved with an 
accuracy, required in modern astrophysical investigations (i.e. at 1 % level) of the recombination history of hydrogen in the early 
Universe. 

We formulated the "two-photon" approximation which make it possible to separate out cascade emission with two-photon 
links and show that this type of cascades gives the contribution to the two-photon transitions, i.e. to the radiation escape, 
comparable with the contribution of the "direct" two-photon transitions. On a basis of this approximation the decay of 3p level 
was described. It was shown that the cascade two-photon decay rates are comparable with the "pure" two-photon contribution 
of the 3s — > Is + 2^{E\) channel. Unlike the 3s — »■ Is + 2"f(El) case, where the contribution of the "pure" two-photon 
decay is nonseparable from the cascade contribution, for 3p — > Is + 37(£1) decay only the cascade contributions with the 
two-photon links should be taken into account. The reason is that the "pure" 3-photon contribution is beyond the accuracy of 
the "two-photon" approximation, adopted in this paper. The four-photon emission process of the 4s level was considered also. 
The result is: the decay rates for the cascades with the two-photon links are comparable with the contribution of the "direct" 
two-photon transitions. 

The main goal of our paper is the formulation of the "two-photon" approximation which allows for the rigarous incorporation 
of all types of the two-photon processes. This may be important for the more accurate astrophysical investigations of the cosmical 
radiation background. 
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2p 

r<Ei> 
— i ls 

Fig. 1. Schematic picture for the transition 2p —¥ ls in hydrogen atom. Horizontal lines denote the electron levels, the vertical line with an 
arrow denotes the photon transition. Transition rate (in r.u.) is (2p — ls) = 732.722ma(aZ) 4 . For hydrogen Z = 1. Here I7 denotes 

the number of photons and El describes the type of the photon. 

ls 
2P 

Fig. 2. The Feynman graph, corresponding to the decay process Fig. 1. The solid line denotes the electron, the upper and lower parts of this 
line correspond to the final and initial electron states. The wave line with an arrow at the end denotes the emitted photon with momentum k, 

frequency oj = |fe| and the polarization e. 

A 

n 

A 



Fig. 3. The Feynman graph describing the elactic photon scattering on an atomic electron. The indices A. n correspond to the initial (final) 

and intermediate electron states. Notations are the same sa in Fig. 2. 
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A 

n l 




n 2 

A 

Fig. 4. The Feynman graph corresponding to the lowest-order electron self-energy insertion into the electron propagator in Fig. 3. The 

internal wavy line denotes the photon propagator in the Feynman gauge. 




A 

Fig. 5. The vacuum polarization insertion in the electron propagator in Fig. 3. The notations are the same as in Figs. 3 and 4. 



2s 



Is 



Fig. 6. Schematic picture for the transition 2s — > Is + 2j(El). Double verticale line with arrow denotes the total two-photon transition (the 

"pure" two-photon) in rhis case. 



A' 



Fig. 7. The Feynman graph for the two-photon emission process corresponding to the A — > A' + 2j transition. All the notations are the same 
as in Fig. 2. Two graphs occur due to the permutation symmetry of the emitted photons. 
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15 



Fig. 8. Schematic picture for the transition 3s — > Is + 2^ (El). Double verticale line with arrow in Fig. 8 a) denotes the total two-photon 
transition, the two ordinary vertical lines in Fig. 8 b) sorrespond to the "pure" two-photon transition and the ordinary vertical lines with 

arrows in Fig. 8 c) describe the cascade photons. The horizontal dashed line in Fig. 8 c) denotes the intermediate energy level. One has to 
remember that Fig. 8 presents the decomposition of the amplitude, so the interference terms between channels 8 b) and 8 c) arise in the 
probability expression. Moreover, as it was explained in the text, actually the contributions 8 b) and 8 c) are inseparable. 



0.025 




Fig. 9. The frequency distribution dWj^ a /dui for the total two-photon transition 3s — > Is + 2-y including cascade and "pure" two-photon 

transitions as functions of the frequency (in a.u.). The values dW^^/duj divided by a 6 (a is the fine structure constant) is plotted versus the 
frequency within the interval [0, uq]. The boundaries for the frequency intervals I-V are also indicated as vertical lines. 
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VIED 



d) 



2p 



Is 



Fig. 10. Schematic picture for the transition 3p — ¥ Is + 3*y(El). Triple vertical line with arrow in Fig. 10 a) denotes the total three-photon 
contribution, the other notations are the same as in Fig. 6. In the "two-photon" approximation only the channels Fig. 10 c) and Fig. 10 d) 

contribute to the decay probability at the adopted level of accuracy. 



A' 



Fig. 11. The Feynman graph for the three-photon emission process corresponding to the A — > A' + 3j transition. All the notations are the 
same as in Fig. 2. The graphs with all permutations of the photon lines should be added. 
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Fig. 12. Schematic picture of the two-photon decay of the state 4s. Notations are the same as in Fig. 7 
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Fig. 13. Schematic picture of the four-photon decay of the state 4s. Quadruple vertical line denotes the total 4-photon contribution. In the 

two-photon approximation only c), d) and e) schemes contribute. 



